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We present the energetic aspect of open systems which may exchange particles with their 
environments. Our attention shall be paid to the scale that the motion of the particles is 
described by the classical Langevin dynamics. Along a particular realization of the stochastic 
process, we study the energy transfer into the open system from the environments. We are able 
to clarify how much energy each particle carries when it enters or leaves the system. On the 
other hand, the chemical potential should be considered as the concept in macro scale, which is 
relevant to the free energy potential with respect to the number of particles. 
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§1. Introduction 

Open systems may exchange matter with their en- 
vironments, in contact with a heat reservoir. The en- 
ergy transfer into such a system from its environments 
may also take place. From the viewpoint of energet- 
ics, the question may arise as to how much energy the 
system gains when a particle enters or leaves the sys- 
tem. In particular, we may ask the question of how 
much energy a particle carries upon migration. For 
instance, we may imagine a macromolecule which can 
bind some small molecule, and may study the varia- 
tion in the energy of the macromolecule upon binding. 1 ) 
According to the thermodynamic relation of open sys- 
tems: AE = W + TAS + fiAn, one may expect that 
each particle which enters into the system carries the 
energy equal to the chemical potential fi. However, we 
should notice that the thermodynamic relation only indi- 
cates the change in the thermodynamic state variables in 
a quasi-static process. In other words, we can evidently 
not discuss the energetic aspect of a particular event of 
particle exchange, based on the thermodynamic relation. 
Indeed, we will show that such an argument is not valid. 
In this Paper, the question will be approached from the 
viewpoint of the microscopic dynamics which describes 
the motion of particles entering or leaving the system. 
From such a viewpoint, we will study the energy balance 
of open systems upon exchanging a particle. We shall see 
that the chemical potential is obtained as the concepts 
in the macroscopic level that we find by considering the 
ensemble of the dynamics. 

In order to study the energy transfer accompanied by 
the migration of the particles, we need a model of the 
system which can express the dynamics of particles 2-4 ) 
and we need to formulate an energetic interpretation of 
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the model. In the present Paper, first we suppose that 
a closed system in contact with a heat reservoir consists 
of N particles. Then, we suppose a subsystem of the 
closed system. The subsystem can exchange the parti- 
cles with the rest of the closed system. Thus, we can 
consider the subsystem as an open system. By studying 
the energy transfer of the open system, we construct the 
energetics associated with the particle exchange. 

Let us first suppose a closed system denoted by fin 
which consists of the N particles. The dynamics of 
these particles is assumed to be described properly by 
the Langevin dynamics (here we only consider the over- 
damped case), given by 

dXi(t) _ dU{x(t),a{t)) 



7- 



dt 



dxi 



•&(*), 



(1.1) 



where Xi is the position of the i-th particle, U(x, a) is 
a potential energy with x = (xi, ■ ■ ■ , xn) and a control- 
lable parameter a, 7 is a friction constant, and £i(i) is 
the white Gaussian processes which is characterized by 
(&(t)) = 0, and (&(ti)&(t 2 )> = 2~ff3- 1 6 itj 6(ti - 1 2 ) with 
= ksT. krj is the Boltzmann constant and T is 
temperature of the heat reservoir. Throughout this Pa- 
per, we consider the processes in which the parameter a 
is supposed to be changed during the time interval from 
t = to t = t by some external system according to a 
given protocol a(t) from a(0) = a,i to a(r) = a/. 

For the energetic interpretation of the Langevin dy- 
namics, we adopt the stochastic energetics. 5 ^ Let 
us consider both the friction and fluctuation terms 
in Eq.(l.l) as the force which the heat reservoir exerts 
on the system. Then, the work done on the system by 
the heat reservoir is given by 
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Fig. 1. Schematic diagram of the present situation. The open system Q consists of the particles being in the spatial domain T> C X>o- 
Hence, Q may exchange particles with its surrounding. How much energy do the particles carry when these enter or leave f2? 



where we use Eq.(l.l) and dU = f£da+ |§dcc. The left 
hand side defines the heat Q[(x(0) 7 ai)^(a;(r), a/)] along 
a particular trajectory. The notation in the bracket [•] in- 
dicates that the quantity takes a particular value which 
is obtained along a particular trajectory x(t) and the 
protocol a(t). We should notice that the above stochas- 
tic integral is interpreted as the Stratonovich integral. 8 ) 
Throughout this Paper, the stochastic integral is used in 
the Stratonovich sense. 

The first term on the right hand side of Eq.(1.2) is the 
change in the internal energy E(t) = U(x(t),a(t)), which 
we denote by AE[(x(0), dj)— >(x(t), a/)]. The second 
term is interpreted as the work TU[(£e(0), ai)^(x(r), a/)] 
done on the system by some external system with chang- 
ing the parameter a(t) from a(0) = ai to a(r) = a/. In 
this way, along one trajectory, the energy balance rela- 
tion 5 ) is obtained as 

A£;[(x(0),a i )-(aj(r),o / )] 

= W[(x(0),a l )^(x(r),a f )} 

+Q[(x(0), ai )^(x(T),a f )]. (1.3) 

This Paper has been organized as follows. A model 
of the open system will be introduced in Section 2. The 
energy balance of open systems along a particular tra- 
jectory is shown in Eq.(2.17). Then we find how much 
energy each particle carries at the boundary of the open 
system. In Section 3, it shall be shown that our study 
of open system is consistent with thermodynamics. It 
will be clarified that the chemical potential is relevant 
to the concept of free energy. Then the amount of work 
onto the open system in the quasi-static process is dis- 
cussed and the thermodynamic relation is obtained. In 
Section 4, we will compare the energy balance relation 
with the thermodynamic relation in Eq.(4.1). Then, wc 
conclude that the energy which a particle carries upon 
migration is not equal to the chemical potential /i. We 
also show some remarks and future problems. 



§2. The Energy Balance Relation of Open Sys- 
tems 

Let us suppose that the particles of the closed sys- 
tem fi are confined in a spatial domain T> , and the 
open system denoted by consists of the particles being 
in the spatial domain D C T>o- Using some function f{x), 
the spatial domain T> is specified by f(xi) < c. Hence, il 
may exchange particles with its surroundings. See Fig.l 
for a schematic illustration of the situation. 

We should note the possibility that no particle is in 
V. In such a case, Q contains no degrees of freedom. 
In order to avoid such a situation, it may be easy to 
introduce another degrees of freedom into f2o, that is, 
a "transducer" which interacts with the particles being 
in T>. Then we can define Q which consists of the trans- 
ducer and the particles interacting with the transducer. 
However, it would not provide further physical insight. 
In order to present our idea simply, here we allow the 
situation that has no degrees of freedom. 

In the following Subsections, we study the energy bal- 
ance relation of the open systems. 

2.1 Open systems as a subsystem of a closed system: 
One particle case 

Before considering many particle systems, we present 
the one particle system, i.e., N = 1. 

Using the Heviside function 0(x), 

««={» <*<»> ■ <"> 

we define the internal energy E(t) of £1 at time t, given 

by 

E{t)=0(c-f{x 1 (t)))U{x 1 {t),a(t)) ( 2 - 2 ) 

We suppose U{x\,a) depends on a if x\ G T>, i.e., 
dU(x\,a)/da = if x\ ^ V. Then, the change in the 
internal energy during the process is given by 

AE[{x 1 (0),a i )^(x 1 {T),a f ))= J dE(t) 
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e(c-f( Xl (t))) 



dU( Xl (t),a(t)) 



dxi (t) 



[ 

Jo 



d9(c-f( Xl (t)j) 
dxi 



U(xi(t),a(t))dxi(t)(2.3) 



The first term on the right hand side is the mechanical 
work VK[(xi(0), ai)^(xi(r), o/)] on done by the exter- 
nal system. The second term is rewritten by substituting 
the Langevin equation Eq.(l.l) as 



0(c -7 



dxx{t) 



+ £i(t) )dx 1 {t). (2.4) 



10 v dt 

Hence, this term is the heat Q[(xi(Q), ai) — > (xi(r),a/)] 
transfered into f2 from the heat reservoir. The third term 
on the right hand side of Eq.(2.3) is transformed into 

- [ T ^^-S(f( Xl ) - c)U( Xl ,a(t))d Xl (t) 

JO Q-X-L 

- df( Xl (t)) y, S(t - t^g) 



dt 



dt 



U(xi(t),a(t))dt 



= U( Xl (t),a(t))- ]T U( Xl (t),a(t)), (2.5) 
t=t in t=t° ut 

where t m and t out indicate all the instants that satisfy 
f{ Xl (t m )) = c,df( Xl (t m ))/dt < and f( Xl (t out )) = 
c,df{x 1 (t out ))/dt > respectively. That is, +U(x 1 ,a) 
is added when f( Xl (t)) — c varies from positive to neg- 
ative and — U{xi 1 a) is added when f{x\{t)) — c varies 
from negative to positive. Hence, the value of this 
term increase or decrease by U(xi,a) when the orbit 
runs across the boundary of V at x\. We call this 
term Q mi a [ ( Xl (0) , a t ) -» (xi (r ) , a f )] . 
Therefore, Eq.(2.3) is rewritten as 

Af?[(a;i(0),Oi) -» (an (t),o/)] 

= W[(x 1 (0),a i )^(x 1 (T),a f )] 

+Q[{x 1 (0),a i )^(x 1 {T),a f )} 

+Q ml °[(x 1 (0),a l )^(x 1 (T),a f )}, (2.6) 

which displays the energy balance relation of the system. 

2.2 Many particle systems with interactions 

Let us next study the many particle systems, i.e., the 

closed system fto contains N interacting particles. 

The interaction potential of the particles is the sum of 

a n-body potential for n — 1, 2, • • •, given by, 

N / N s 

u( X ,a) = j:(€ ) (-)+ E OS, + ■•■)) ( 2 - 7 ) 

ii = l ^ J2=ii+1 ' 



where 



potential of , ■ ■ ■ , iEi n . We assume that the one body 
energy potential (jfp(a) = <fi^ (xi, a) depends on a if 
Xi e V, i.e., d<j)^(a)/da = if Xi £ V. We assume only 
short range interactions so that the n body potential has 
non-zero value only when all the n particles are close to 
each other. 

As in the case of the one particle system, the internal 
energy E of ft is defined by considering whether each 



particle is in V or not: 

N 



N s 

+ E -))). 

(n) 

where 9^ ... i is a function indicating whether the n par- 
ticles belong to ft or not. Using the Hcvisidc func- 
tion 9(x), 9\™\.. in = 9 { - n \x il , ■ ■ ■ ,x in ) is defined by 

Cl^i-a-C^a-C)' ( 2 -9) 

The value of 9 i± ... in is unity when at least one of the 
n particles is in Q, or is zero, otherwise. This means 
that certain n particles belong to when at least one 
particle of these particles is inside V. Thus, 9^\.. in = 1 

and the interaction energy (j>^\.. i of the n particles is 
included in the internal energy E. However, if m(< n) 
particles ji, • • • , j m of these n particles are outside V, 
9j™J. j m = and the interaction energy ( t ) j\\.j m of 
these m particles is excluded from the internal energy E. 

The change in the internal energy along a particular 
trajectory is given by, 



AE[(x(0), ai )^(x(r),a f )} 



V),^ 1 ' («(*)) 



= E/ *< (1) w 
+E / 



9a 



da(i) 



-7^+&(t))<W0 



, jo 



o 9x 



E 



(2.10) 



with 



N N 

= # } +E(«+E(C 



+ 



ji=i 



J2=jl + 1 

i2 5^« 



+ E 



lJl,-'JIV-l 



(2.11) 



</> (n) (ajj! , • • • , x in ) is the n-body energy an d 



JN-l=jiV-2 + l 



ii=i 

N 



J2=jl + 1 
J2^» 



e ei.cj - (2-i2) 



JiV-l=JJV-2 + l 

JJV - 1 #i 
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Eq.(2.10) is obtained by the following transformation. 
The total derivative of the internal energy Eq.(2.8) is 
written as 

dE = ye^^M d a 

1 da 

i—l 

N A / N N 

E(4 , #+E(&&+- 

2=1 1 jl=l 32=31 + 1 

JN=jN-l + l 

Considering that ft,-^ „■ and (/i^ „■ are symmetric 
functions of Xj x , • • • , Xi N , the second term on the right 
hand side of Eq.(2.13) is rewritten as 

N Q / N N 

+ U e SM + E (CL^SL+- • • 

1=1 * V Jl=l J 2 =il+1 

•••+ E C.^-.C,-,^)))^ (2-14) 

iw-l=jiV-2 + l 

Noting that 0,- n ^ , is transformed as 0^ ,• = tf,- 1 ^ + 
(1 - 6 , - 1 1) )6'-"7 1 j n , then Eq.(2.14) is rewritten as 

i—i 

with Eqs.(2.11) and (2.12). Using the relation = 
-7^- + Eq.(2.15) is rewritten as 

«— i x / 

Substituting Eq.(2.16) into Eq.(2.13), and integrating it 
with respect to x(t) and a(t) from t = to t = r, then a 
simple calculation leads to Eq.(2.10). 

The first term on the right hand side of Eq.(2.10) is the 
work VF[(a:(0), ai)-^(x(r), a/)] on fi done by the external 
system. The second term of Eq.(2.10) is the heat trans- 
fered from the heat reservoir to the particles being in T>, 
denote by Q[(x(0), ai)— >(x(t), ffl/)]. As is discussed in 
the previous Section, the third term of Eq.(2.10) equals 
to <pi{t m ) when the orbit Xi(t) enters into V over the 
boundary at t = t m , or equals to —(pi(t out ) when the 
orbit Xi(t) goes out of V at t = t out . Let us denote 
this energy transfer accompanied by the migration of 
the particles by Q mi a[(x{0), (h)->(x(t), a f )]. The fourth 
term of Eq.(2.10) distinguishes the interacting particles 
systems from one particle systems or free particles sys- 
tems. This term shows the contribution by the particles 
outside T> which interacts with the particles inside V 
across the boundary dV of T>. Here, we denote this term 
as Q dv [{x(0),ai)^(x{T),a f )]. 

Therefore, Eq.(2.10) is rewritten as 

AS[(x(0),o i )-.(a : (r),o / )] 
= W[(x(0), ai )^(x(T),a f )} 



+Q[(x(0),ai)^(x(r),a f )} 

+Q mig [(x(0), ai )^(x(r),a f )} 

+Q dv [(x(0),a^(x(T),a f )] (2.17) 

In this way, along a particular trajectory, we have the 
energy balance relation of the open systems with many 
body interactions. 

In the similar way, the change in the ensemble average 
of the internal energy E{t) can be also discussed. The 
ensemble average of the internal energy is defined by 

(E(t)} = J E(t)P(x,t)dx, (2.18) 

where P(x,t) is the probability distribution function of 
the state of fio obtained by solving the Fokker-Plank 
equation, and (•) indicates the ensemble average over all 
the possible trajectory. 

The energy balance relation of the mean values is 

A(E[ ai ^af}) = (W[a t ^a f }) + (Q[a z ^a f }) 

+ (Q m9 [a l ^a f }) + {Q dv [a t ^a f ]). (2.19) 

Henceforth we use the notation [a%— >af\ for those quan- 
tities concerning the processes in which a(t) is changed 
from di—tCLf. Each term on the right hand side corre- 
sponds to each term in Eq.(2.17). We do not need to 
postulate further condition in order to obtain this ex- 
pression for the mean values. Comparing Eq.(2.17) with 
the thermodynamic relation that we have presented in 
Section 1, we have 

(Q) + (Q mi9 ) + (Q dV ) = TAS + nAn, (2.20) 

where we suppose the quasi-static process. We will dis- 
cuss the meaning of this expression in Section 4. 

Before concluding this Section, let us evaluate Q mtg 
for the case that we are able to effectively exclude the 
possibility of the direct interaction across the bound- 
ary dV, i.e. Q 9V — 0. 9 ) We further suppose U(x, a) — 
U(x,a) defined by 

U(x,a) = U {n) { Xl , ■ ■ ■ ,x n ,a) + (N - n)U , (2.21) 

when Xi satisfies /(a;,) < c, (i = 1, • • • , n) and Xj satisfies 
f(xj)>c, (j = n + 1, • • • , N). This assumption leads to 
the simple form of Q mig as 

Q mig = ^o(E 1 - E x )- ( 2 - 22 ) 

t=t in t=t° ut 

Even in the present simple case, the relation between the 
terms on the left and the right hand side of Eq. (2.20) does 
not seem evident. However, this does not imply that our 
study of the open system based on the stochastic ener- 
getics is zrzcompatible with thermodynamics. In the next 
Section, we will show the consistency. Then in Section 4, 
Eq.(2.20) is consistently interpreted mentioning the dif- 
ference of the scale of description on each term. 

§3. The thermodynamic relation from micro- 
scopic dynamics 

It has been shown in the previous works 5-7 -* that the 
framework of the stochastic energetics is consistent with 
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thermodynamics of the system in contact with a heat 
reservoir. In this Section, we shall show that even in the 
open system the stochastic energetics is consistent with 
thermodynamics. 10 ) First, in Subsection 3.1 we study 
the probability that the open system f2 is found to be a 
particular physical state which is specified by the num- 
ber of particles in T> with their positions. Then, we will 
find that the chemical potential is relevant to the concept 
of free energy potential obtained through the reduction 
of the description of the system. In Subsection 3.2, we 
also study the quasi-static process, in which the work 
onto the open system is shown to equal the difference 
of a pertinent thermodynamic potential. Then, in Sub- 
section 3.3 the thermodynamic relation is obtained so 
that the compatibility of our study with thermodynam- 
ics shall be shown. 

In the present Section, the state of the closed sys- 
tem f2 is specified by x^ N > = (x\, ■ ■ ■ ,Xn) which con- 
stitute the phase space T^ N \ We assume that all the 
particles are identical. If two points in differ only 
by the permutation of x^s they may be considered as 
representing the same physical state. Then, we can intro- 
duce the phase space r'W obtained by the identification 
of all such points. In r'W the state of the particles is 
specified by x'( N \ For each point x'W in T'( N \ there 
are N\ points in T^ N \ Hence, if f(x^) is a symmetric 
function of all the Xj's then 

/ f( X >W)dx'W = ±- [ /(*W)d*W 
J(v N )' AM Jt> " 

(3.1) 

where (Dq N )' indicates the domain in r'W which cor- 
responds to V a N in rw. 

3. 1 The chemical potential as free energy potential 

The equilibrium distribution function of of Q is 
obtained by solving the Fokker-Plank equation, 



Hence, U^^x'^) is given by 



P eq {x^ N \a) = 



1 1 



Z{a) v» 

with the normalization condition: 

1 1 f airl^im 



-PU(x( N \a) 



1 = 



Vq Z{a) J (t>0 n 



e -PU(X^\a) dx {N)_ 



(3.2) 



(3.3) 



Here the constant vq of the dimension of volume has 
been introduced so that Z{a) is dimensionless. Since 
U{x^ N \a) is symmetric function of x^s, the probability 
distribution function of the state x'W is 

with sb'W e {vgy. 

By integrating P' eq (x 1 ^ , a) with respect to the de- 
grees of freedom of the particles being outside V, that 
is, in T>o\D, we obtain the probability distribution func- 
tion of that there are n particles in V with their posi- 
tions jc'M G {V n )' as 

1 1 -pu <n) (X' <n) ,a) 



Z{a) < 



(3.5) 



where U^ n \x'^) is considered as the free energy po- 
tential under the restriction that the state of fl is x'( n \ 



Nl 

N-r 



e -0U{X'^>,a) dx l{N-n)^ 

((T> \T>)»-ny 

(3.6) 

with G (V n )'. We can show the following normal- 
ization condition holds: 

JV 

1 = E/ P& n) {x Kn) ,a)dx'W. (3.7) 

Hereafter we shall exclude the possibility of the direct 
interaction among the particles across the boundary dV 
and assume U(x^ N \a) = U{x^ N \a), with Eq.(2.21). 
Then, Eq.(3.6) is rewritten as 

e -0U^\x'^\a) 

= N\(V - V) N ~ n (uin Hx K nKa)+(N _ n)Uo ) 

v£- n (N-n)\ ' V ' ' 

with jc'(") G (£>")'. We shall consider the case that the 
volume of Qo ; Vq = Jt> ^ x an ^ N are sufficiently large 
with p = N/Vo and V = J v ldx being kept to be finite 
constants. Then, U^ n \x'^ n \a) is given by 

W W (x„,a) = U^\x'^\a) - n(U + log pv ) 

+N(u + (3- 1 log^) +/3- 1 pV + o(±). (3.9) 

Since, on the right hand side, only the first and the sec- 
ond terms depend on x'( n > and on n, these terms should 
be regarded as a new free energy potential assigned to the 
situation where there arc n particles in T> with their po- 
sitions corresponding to the physical state x'^ G (D n )' . 
Hence, introducing the chemical potential p through the 
definition: 



p = U + P 1 \ogpv , 



(3.10) 



Peq 1 ^ {x 1 ^ , a) is rewritten as 
P#V< n >,a) = 1 i. e -/5(^ ) (x'W,a)-^) ; (3 n) 

where S(a) is a normalization constant given by 
3(a) = lim Z(a)e fm ( U ^ llo ^ + °^) 

oo . - 

Eq.(3.11) is the grand canonical distribution. 

Notice that the chemical potential p is obtained as 
a part of the free energy potential U^ n \xi, ■ ■ ■ ,x n ,a). 
Hence, the crucial point for obtaining the chemical po- 
tential is to consider the trajectory x'( N ^ in the reduced 
space r'W rather than and to extract the effective 
potential of the state x' < - n ' G (T> n )' by eliminating the 
degrees of freedom in ((Do\D) N ~ n )' . 

3.2 Work on the open system in quasi-static process 

As is discussed in the previous Section, the work onto 
the open systems is a random variable, whose value 
is determined depending on a particular trajectory. If 
the process is quasi-static, however, we expect that the 
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amount of work is deterministic quantity. In this Sub- 
section, we show that the work onto the open system 
is given by the difference of a pertinent thermodynamic 
potential. 

Here we note that the Langcvin dynamics Eq.(l.l) has 
the ergodic property that long time averages can be re- 
placed by the ensemble averages. When the change of the 
parameter a is slow enough and hence |da(i)/di| is small 
enough, the work along a particular trajectory given by 
f£ du ^ x d ^' a ^ da(t) can be evaluated asymptotically 
as the integral of the equilibrium ensemble average of 
the integrand with respect to a from aj to a/. Using 
P' e ^ (x'( n ) , a) , the work done by the external system is 
written in the limit r^ooas 

lim W[(x(0), ai )^(x(T),a f )] 

T^OO 

= AJ[oj-»o/] = J{a f )- J(ai), (3.13) 

where J (a) is the grand canonical potential defined by 
J(a) = -/3" 1 logS(a). 

Hence, if the change of the parameter is slow enough, 
the work done by the external system onto the open sys- 
tem equals to the difference of the grand canonical poten- 
tial with the probability being unity, without mentioning 
the ensemble average. 

It should be noted that for a finite time interval, the 
amount of the work of a particular process can be larger 
or smaller than AJ. Considering the ensemble average 
of the work, the minimum work principle holds on the 
open system, 10 ^ i.e., (W[aj— >a/]) > AJ[aj— »a/]. 

3.3 Thermodynamic relation 

The equilibrium ensemble average of the internal en- 
ergy of fl is given by 

N 

(E(a)) = Y / U^\x'^)P^\x'^\a)dx'^ 

= E / (-0~ 1] °S (vZPfrHx'W, a)) + J(a)) 

xP^(x' (n \a)dx' (n \ (3.14) 

where we substitute Eq.(3.11) into U {n \x l( - n \a). Intro- 
ducing the entropy S{a) given by 

N 

S(a) = -k B J2 P& n) (vZP>M)dx'M, (3.15) 

n=0- / ( CB )' 

we have 

(E(a)) =J(a)+TS{a) + fj,(n{a)}, (3.16) 

where (n(a)) is the equilibrium ensemble average of the 
number of particles in V. 

When the parameter a is changed from a» to ay quasi- 
statically, the change in the internal energy is given by, 

A(E[ai -> a/]) = W[a,i -» a f ] 



+TAS[ ai -> a/] + ^A(n[a; -» a/]), (3.17) 

where Eq.(3.13) is substituted. In this way, we have the 
thermodynamic relation of the open system, which holds 
on the ensemble of the trajectories x(t). 

Eq.(3.17) may suggest that the increase in number of 
particles during a quasi-static process would accompany 
the transfer of energy into the system as if each particle 
carries the chemical potential. In the next Section, we 
conclude this Paper, mentioning how the thermodynamic 
relation is compatible with the energy balance relation 
that we have obtained in Eq.(2.17). 

§4. The energy carried by a particle upon mi- 
gration: Conclusion 

In the preceding sections, we have constructed the en- 
ergy balance relation Eq.(2.17) of the open system focus- 
ing our attention to the scale that the system is described 
by the Langevin dynamics. On the other hand, it has 
been clarified that the chemical potential, which reflects 
the probability that the open system contains a partic- 
ular number of particles, is relevant to the free energy 
potential. Here, we shall illustrate the compatibility of 
the energy balance relation Eq.(2.17) with the thermo- 
dynamic relation Eq.(3.17), noting the difference of the 
scale of their description. 

As is the case in the previous Section, here we shall 
consider the case that the possibility of the direct in- 
teraction among the particles across dT> is effectively 
excluded, and we suppose the quasi-static limit and 
V, N^oo. From Eqs.(2.17) and (3.17) we have the quan- 
titative relation 

(Q) + (Q mi9 )=TAS + ^A{n), (4.1) 

and Q dv = 0. If we assume U(x,a) = U(x,a) with 
Eq.(2.21), we have Q mi <> = U An. Remembering the def- 
inition of fj, given by Eq.(3.10), we have, 

<Q"^)^ M A<n>. (4.2) 

This leads to our main conclusion of the present Paper 
that each particle does not carry jj, across the boundary 
of the system, when it enters or leaves the system. 

We should notice the difference between the energy 
balance relation associated with the left hand side of 
Eq.(4.1) and the thermodynamic relation associated with 
the right hand side of Eq.(4.1). The thermodynamic re- 
lation is a closed expression in terms of the variation in 
the thermodynamic state variables of a system. On the 
other hand, the energy balance relation is the relation 
of energy exchange with the surroundings. Let us study 
how (Q) is expressed by using the thermodynamic state 
variables. From Eq.(4.1) and Q mig = U An, we have 

(Q) = TAS* + k B Tlog (pv )A(n) ? TAS. (4.3) 

This quantitative relation shows that the ensemble av- 
erage of the energy transfer from the heat bath is ex- 
pressed by the variation in the two thermodynamic state 
variables AS and A(n) in both the two terms TAS and 
juA(n) of the thermodynamic relation. 

We should also mention the difference of the level of 
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the description between the energy balance relation and 
the thermodynamics relation. Since the energy balance 
relation holds on a particular trajectory, it is not neces- 
sary for obtaining the relation to know the probability 
of the realization of the trajectory. In this sense, the 
energy balance relation describes the microscopic scale. 
On the other hand, in order to obtain the thermody- 
namic relation, it is necessary to know the probability 
distribution of the physical states. The thermodynamic 
relation is obtained so that it is a relation among the 
ensemble averages. In this sense, the thermodynamic 
relation describes the macroscopic level. Consequently, 
the chemical potential is regarded as the concept in the 
macroscopic description. 

Thus, when we observe a particular process, i.e. a 
sample path, we can not find the chemical potential as 
well as the thermodynamic relation, but measure a par- 
ticular value of Q obeying the energy balance relation. 
With the increase of the number of the samples, many 
different values of Q shall be recorded. Then, we know 
that Q takes a value following some probability distri- 
bution. Eventually, we find the chemical potential and 
the thermodynamic relation by examining the ensemble 
of the processes. 

It may be a future problem to study an energy conver- 
sion cycle of a micro-machine with at least two particle 
reservoirs. For example, suppose that the two reservoirs 
consist of the particles of the same species and provide 
the same one-body energy potential and the same tem- 
perature environment. Then the difference of the chem- 
ical potentials fin and fi^ of these reservoirs can arise 
only through the difference of the density of the parti- 
cles. Thermodynamics tells that the maximally avail- 
able mechanical work obtained per particle transported 
is /iff — hl (> O). 4, n > By studying this process with the 
micro-machine, it might be more evident how the small 
fluctuating system feels the chemical potential so as to 
produce the mechanical work. 
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